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Abstract 

The approximation of matrices to the sum of tensor products 
of Hermitian matrices is studied. A minimum decomposition of 
matrices on tensor space Hi ® H2 in terms of the sum of tensor 
products of Hermitian matrices on Hi and H2 is presented. From 
this construction the separabihty of quantum states is discussed. 
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1 Introduction 



The quantum entangled states have become one of the key resources in quantum information 
processing. The study of quantum teleportation, quantum cryptography, quantum dense 
coding, quantum error correction and parallel computation [H |21 E] has spurred a flurry of 
activities in the investigation of quantum entanglements. Despite the potential applications 
of quantum entangled states, the theory of quantum entanglement itself is far from being 
satisfled. The separability for bipartite and multipartite quantum mixed states is one of the 
important problems in quantum entanglement. 

Let Hi (resp. H2) be an m (resp. n)-dimensional complex Hilbert space, with 
i = 1, m (resp. j = 1, n), as an orthonormal basis. A bipartite mixed state is said 
to be separable if the density matrix can be written as 

P = ^Pipl ® p1, (1) 
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where < j9j < 1, YliPi ~ Pi ^^"^ pf rank one density matrices on Hi and H2 
respectively. It is a challenge to find a decomposition like ((T)) or proving that it does 
not exist for a generic mixed state p |5 El E] . With considerable effort in analyzing the 
separability, there have been some (necessary) criteria for separability in recent years, for 
instance, the Bell inequalities [Zj, PPT (positive partial transposition) jHI (which is also 
sufficient for the cases 2x2 and 2x3 bipartite systems [Oj), reduction criterion [TIH ITT], 
majorization criterion [I^j, entanglement witnesses jH] and realignment [15^ it^ IT7j 

and generalized realignment |TH|, as well as some necessary and sufficient criteria for low 
rank density matrices [T9 | 120 } l2T] and also for general ones but not operational [Oj. 

In |221 the minimum distance (in the sense of matrix norm) between a given matrix and 
some other matrices with certain rank is studied. In [23] and j21], for a given matrix A, the 
minimum of the Frobenius norm like ||^ — -Bi ® Cj||ir is investigated. In this paper we 
develop the method of Hermitian tensor product approximation for general complex matrix 
A, i.e. we require Bi and Ci to be Hermitian matrices. By dealing with the Hermitian 
condition as higher dimensional real constraints, an explicit construction of general matrices 
on Hi®H2 according to the sum of the tensor products of Hermitian matrices as well as real 
symmetric matrices on Hi ® H2 is presented. The results are generalized to the multipartite 
case. The separability problem is discussed in terms of these tensor product expressions. 



2 Tensor product decomposition in terms of real sym- 
metric matrices 

We first consider the tensor product decompositions according to real symmetric matrices. 
Let A be a given mnxmn real matrix on Hi®H2. We consider the problem of approximation 
of A such that the Frobenius norm 

r 

\\A-Y,B^®C,\\p (2) 

i 

is minimized for some m x m real sjTiimetric matrix Bi on Hi and n x n real symmetric 
matrix Ci on H2, i = 1, r G IN. 

We first introduce some notations. For an m x m block matrix Z with each block Zij of 
size n X n, i, j = 1, ...,m, the realigned matrix Z is defined by 

Z = [vec{Zii), ■ ■ ■ , vec{Zmi), ■ ■ ■ , vec{Zim), ■ ■ ■ , vec{Zmm)f, 
where for any m x n matrix T with entries tjj, vec{T) is defined to be 

Vec{T) = [til, ■ ■ ■ , tml,tl2, ■ ■ ■ , tm2, " " " , ^In, " " " , ^mn]*- 

There is also another useful definition of Z, {Z)ij^ki = {Z)ikji. A matrix Z can be expressed 
as the tensor product of two matrices Vi on Hi and V2 on H2, Z = Vi ® V2 if and only if 
(cf, e.g., I21|) Z = vec(yi)vec(y2y , i.e., the rank of Z is one, r{Z) = 1. 

Due to the property of the Frobenius norm, we have 

r r 

\\A-Y,Bi(^Ci\\F=\\A-J2^(^<B,)vec{Qy\\F. (3) 

i=l i=l 
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The symmetric condition of the matrices Bi and Cj can be expressed in terms of some real 
matrices 5*1 and 5*2 in a form 



Slvec{Bi) = S'^veciCi) = 0, 



1 T 



(4) 



We define Qs to be an x matrix such that, if we arrange the row indices of 

Qs as {11, 21, 31, ml, 12, 22, 32, m2, mm}, then all the entries of Qs are zero except 
those at 21 and 12 (resp. 31 and 13, ...) which are 1 and —1 respectively in the first (resp. 
second, ...) column. We simply denote 



Qs — [{^21, —612}; {631, —613}; {em,m-l, —Gm-l,m}], 



(5) 



where {e2i, —612} is first column of Qs, with 1 and —1 at the 21 and 12 rows respectively; 
while {631, —613} is second column of Qs, with 1 and —1 at the 31 and 13 rows respectively; 
and so on. 



Similarly we define Qa to be an x 



m(m+l) 



matrix such that 



Qa = [{en}; {621, 612}; {631, 613}; ...; {622}; {632, 623}; {642, 624}; •••; {e^.m-i, em-i,m}, {emm}]- 

(6) 



For m = 2, we have 



Qs 



( \ 

1 

-1 
V ) 



Qa 



/ 1 \ 

1 

1 
Vo 1/ 



Si can then be expressed as, something like QR decomposition. 



Si — Qs = Qi 







(7) 



where Ri is a full rank x matrix, Qi is an orthogonal matrix, Qi = [QsQa) , 

where Qs and Qa are obtained by normalizing the norm of every column vector of Qs and 
Qa to be one. 



For the case m 



Qi 



/ 1 \ 
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V 001/ 



Ri 



V2]. 



S2 has a similar QR decomposition with S2 = Q2 ( ^ ) ; by replacing the dimension 



m with n in ((Tj). 
Set 



Q{AQ2 



An 

A21 A22 



(9) 



Suppose the singular value decomposition of A22 is given by ^122 = J2l=i V\uiv\, where r 
and A,;, i = 1,2, ...,r, are the rank and eigenvalues of ^22^22 respectively and Ui (resp. Vi) 



are the eigenvectors of the matrix ^22^22 (resp. J^^fi-^'^'i) ■ Set Ei = y/XiUi, Ci = Vi. 

[Theorem 1]. Let A be an mn x mn real matrix on Hi ® H2, where dim{Hi) = m, 
dim{H2) = n. The minimum of the Frobenius norm \ \A — Yll -^i ® C'jIIf is obtained for 
some m X m real symmetric matrix Bi on Hi and n x n real symmetric matrix Ci on H2, 
given by 



vec{Bi) = Qi 

[Proof]. Set 

Q\vec{B,) = 
From (0]) and ((Tj) we have 

{Ri 0) 







vec{Ci) = Q2 



Ql vec{Ci) 







Ci 



B,; 
B, 



0, {Ri ) 



0, 



which give rise to Bj = Cj = 0, due to the nonsingularity of R{ and R^- 
From 0, Q and ^ we obtain 

r r 

\\A -Y^Bi(»Ci\\F = \\Q{AQ2 -J2QW<Bi)vec{CiyQ2\\F 

i=l 1=1 

All Ai2\ f Q 



A21 A22 

All A12 
A21 A22 



E 

i=l 

r 

i=l 



B, 



(0 C\) 




i3,,d* 



\A 



+ \\Ai2\\l+\\A2i\\l + \\A22-Y,BA 



t\\2 
i\\F ■ 



(10) 



(11) 



i=l 



From matrix approximation we have that ^22 = X]i=i BiC\ is the singular value decomposi- 
tion (SVD) for A22, which results in formula (fTUI) . □. 

From Theorem 1 we see that if a real symmetric matrix A has a decomposition of tensor 
product of real symmetric matrices, then An = A12 = A21 = 0. As an example we consider 
the Werner state 
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where |^ ) = (|01) — |10))/-\/2. State pw is separable for F < 1/2 and entangled for 
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1 /2 < F < 1 . According to the definition of realignment we liave 
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Here tlie dimension m = n, lience Q2 = Qi is given by (jHI). From we liave 



QiPwQi 
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Pw21 Pw22 
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Tlierefore p^ is generally not decomposable according to real symmetric matrices because 
(Pii')ii = (4-^ ~ l)/6 7^ as long as F 7^ 1/4. From the singular value decomposition of 

(P«>)22, 

\ 

^ _ 6 ' 

(Pu>)22 



2F+1 



V 



6 







6 



1-F 
3 



we have: 



Ml = Vi 




eu2 = V2 




eu-i = V3 




with eigenvalues Ai = 1/4, A2 = A3 = (1 —4F)^/36 respectively, where e = (1 — 4F)/|1— 4F|. 
From (Uni) we have vec{Bi) = ^/X^{l/V2, 0, 0, l/v^)*, vec{B2) = e^/)^{-l/V2, 0, 0, 1/^2)*, 
vec^B^) = e^/\3{0, l/-\/2, 0)*. Therefore the best real symmetric matrix tensor prod- 
uct decomposition is 



1 



Pv 



'2x2 



where are Pauli matrices ai 
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3 Hermitian tensor product decomposition of Hermi- 
tian matrices 

We consider now the tensor product decompositions according to Hermitian matrices. Let 
A be a given mn x mn complex matrix on Hi ® H2. We first consider the problem of 
approximation of A such that the Frobenius norm 



\A-B®C\\f = \\A- vec{B)vec{Cy\ \f. 



(13) 
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is minimized for some mx m Hermitian matrix B on Hi and nx n Hermitian matrix C on 

In order to impose the Hermitian condition of the matrices B and C, we separate the 
matrices B and C into real and imaginary parts such that B = B + iB, C = C + iC, where 
B and C (resp. B and C) are the real (resp. imaginary) parts of B and C respectively. As 

vec{B) = t>ec(B) + ivec{B), we have 



vec{B)vec{Cy 



f ec B fec C 



vec{B)vec{CY) + i{vec{B)vec{CY + vec{B)vec{c) 



We now map the complex matrix A to be a real one: 

A — > 



A A 

-A A 



where A and A are the real and the imaginary parts of A respectively. Now the approx- 
imation problem of complex matrices to the tensor product of two Hermitian matrices is 
reduced to the problem of real matrices and the results in [23 121] can be used accordingly. 
The problem to minimize \\A — vec{B)vec{Cy\\F is reduced to minimize 



A^ A 
-A A 



fec(B) vec{B) 
-vec{B) fec(B) 



fec(c) —vec{C) 
vec{C) fec(c) 



(14) 



under the Hermitian condition: B = B\ C = C^, i.e., B and C are symmetric, B and C are 
antisymmetric. This condition can be expressed in terms of some real matrices 5*1 and 5*2 
in a form 



t ( wec(B) 



Si 



f ec(c) 
±f ec(C) 



0. 



[Lemma 1]. Condition (fTH|) has a QR decomposition such that 



Si = Qi 



Ri 




82 = Q2 



R2 





(15) 



(16) 



where Ri and R2 are full rank matrices, Qi and Q2 are orthogonal matrices. 
[Proof]. Si can be generally expressed as 



^1 



Qs 






Qa 



where Qs and Qa are given by (0) and (jHI) respectively. The QR decomposition of is 
given by 



Qi 



Xi Yi 
Yi Xi 
.2 



(17) 



Q Qa 
Oa 

with Qs and Qa given in section 2, Xi (resp. Yi) is an m? x m? matrix with the first 
m{m — l)/2 (resp. last m{m + l)/2) columns replaced by the matrix Qs (resp. Qa) and the 
rest entries zero, Ri is a diagonal matrix with diagonal elements either 1 or \/2. For the 
case m = 2, 



Xi 



( 

1 

\_ 

/ 



\ 





Yi 



/O 1 





\ 





1 

V2 

1 

V2 








1 / 



Ri 



/ V2 \ 

10 

V2 

\ 1 / 



(18) 
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5*2 has a similar QR decomposition with 

by replacing the dimension m with n in the expression of 5*1 . □ 
Set 

, ( t;ec(B) \ / B \ / t;ec(B) \ / B \ 

(20) 

, / t;ec(c) \ / C \ / t;ec(c) \ C \ 

^2 ^ -vec{C) )-\-c)' ^2 ^ t;ec(C) )-\c)- 

From (P3|) and (fTB|) we have 

{R\ o)(^^)=0' (^i o)(|)=0' 
(i?^ o)(^^^)=o, (i?* o)(^)=o, 

which give rise toB = B = c = C = 0, due to the nonsingularity of R\ and R\. 
Therefore we have 

t;ec(B) \_^( B \_(X, Y,\( \ _ / 
-t;ec(S) ) "^'[-B J [Y^ X, J[-B J [ -Xj 

vec{B) \_^fB\_fX, Y,\fO\_fY,B 
vec{B) J "^'{b J [y^ J{b J {X,B 

Thus -YiB = YiB, XiB = XiB and 

B=iX, + Y,)-\X^ - Y,)B =1 I ( -Qa)B = I^J, (22) 



(21) 



where /™a = diag{I^, — /™), /™ (resp. /™) is an m(m — 1)/2 (resp. m(m + l)/2) dimensional 
identity matrix. 



Let P denote the permutation matrix, 

P = 



I-m'^ X 



It is easily seen that PQiP = Qi- From the second formula in ()21|) we have 



Hence we have 



vec{B) ] ^ ( 13 
vec(B) / ~ I 



t . vec{B) vec{B) \ _ f B , 
^1 ^ -vec{B) vec{B) J ~ [ -B ' ' ^ ' 
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and, similarly, 



where 



Q2 



t f fec(c) —vec{C) 
vec{C) fec(c) 



-C 
C 



(24) 
(25) 



= diag{Ig, —la), (resp. /^) is an ?7,(r;, — l)/2 (resp. n(n + l)/2) dimensional identity 
matrix. 



Set 



That is 



A A 
-A A 



Q2 



All A12 
A21 A22 



(26) 



All = XiAX2 + YI~aY2 + X{AY2 - YIAX2, 
A12 = X{aY2 + YIAX2 + XIAX2 - YIAY2, 
A21 = YIAX2 + X{aY2 + YiAY2 - XIAX2, 
A22 = YIAY2 + X*aX2 + YIAX2 - X{AY2. 

[Theorem 2]. To minimize is equivalent to minimize the following formula 



\M2^BC^\\l^\\Aii^BC^\\l 



\Ai2\\\^\\Ml\\\. 



(27) 



[Proof]. From (dH), ((221) and (j21I) we obtain 



A A 
-A A 



vec{B) vec{B) 
—vec{B) fec(B) 



f ec(c) 
vec{C) 





A 




A 


All 


A12 


-A21 


A22 


All 


A21 


-A21 


A22 



Q2 - Q\ 



B 
-B 

-BC' 




fec(B) vec{B) 

-vec{B) fec(B) 

-C^' 

. C 



-BC' 



—vec{C) 
vecic) 



fec(c) — fec(C) 
vec{C) fec(c) 



\A22 + M'\\l + \\Aii + BC'\\l + \\Ai2\\l+\\A2i\\l. 



□ 

[Lemma 2]. If the matrix A = A + iA is Hermitian, we have the relations: 

A12 = Ml = 0, All = IZ^22i:,a- 

[Proof]. As the matrix A = A + iA is Hermitian, i.e., A is symmetric and A is antisym- 
metric, we have 



{^)ij,kl — {^)ik,jl — {^)il,ik — (yA)ji,lk, {A)ij^kl — iA)ik,jl — —{A)jl^ik — — (^)jj,«fc- 
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Noting that in our construction of and Y^, {Xa)ij^ki 
a = 1, 2, we obtain 



'(^-^a)ji,kh i,^a)ij,kl (Ya^ji^kh 



{X^AY2)ij^pq — {Xi)kl,ij{A)ki 

mn,pq 



0. 



Similarly we have F/AXa = XIAX2 = F/^K, = 0. Hence = ^21 = 0. 

From the relations J™ Fi = -Fi, I^^Xi = Xi and Fa^^a = -^2, 4"a^2 = ^2, we have 
/™ (F/AF2 + X*aX2)JI:, = YIAY2 + X*aX2 and I^,{YlAX2 - X^F^)/!;. = -{YIAX2 - 
XiAY2). Therefore in = ITaA22lsa or A22 = ITa^ul 



n 

s,a' 



□ 



For Hermitian matrix A, by using ()22|) . and Lemma 2 we have \\Aii + -BC*||ir = 
||in + I^^BC^imF = ||/I"„in/"a + ^C*\\f = \ \A22 + BC^Wf- From Lemma 2 we have that 
to minimize | |y4 — i? Cg> C| |f (fTSj) is equivalent to minimize \ \A22 + -BC*| which maybe zero, 
when A22 = -BCK 

Now the minimum of the Frobenius norm \ \A — Xll -^i ® can be obtained readily. 
Suppose the singular value decomposition of A22 is A22 = Y^i=i V^iUivj, where r and A,, 
i = 1,2, ...,r, are the rank and eigenvalues of ^22^22 respectively and Ui (resp. Vi) are the 
eigenvectors of the matrix ^22^22 (resp. 0422^22)- Set Bi = y/Kui, Ci = —Vi. By using the 
results in |23| , for Hermitian matrix A the minimum of 1 1 A — ^^^-^ (g) Cj 1 1 is obtained 
wheni22 = -E:=i^4. 

[Theorem 3]. Let A be an mn x mn Hermitian matrix on Hi ®H2, where dim{Hi) = m, 
dim{H2) = n. The minimum of the Frobenius norm | |A — -Bj ®Cj| |i? is obtained for some 
mxm Hermitian matrix B on Hi and nxn Hermitian matrix C on H2, if A22 = — Yll=i ^i^h 
where ^22 is defined by = Bi + iBi, Cj = Cj + iCi, are given by the relations 



vec{Bi) 
—vec{Bi 







vec{Ci) 
vec{Ci) 



Q2 







:28) 



As an example we consider the bound entangled state on 2 x 4 
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2, n = 4) [23, 
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(29) 



where < 6 < 1. pb is a PPT but entangled state. The QR decomposition in our case only 
depends the dimensions. Qi is still given by (fTTj) and (fTK|) . Q2 is a 32 x 32 matrix with X2, 
F2 in (Uni) given by X2 = (fi, f2, fs, fs, fe, Oio), F2 = (Oe, ai, a2, as, a^, as, as, a^, as, ag, aio). 
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where s,- and a,- are 16 x 1 column vectors: 



fi = (0, 1/^2, 0, 0, -1/^2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) 
f2 = (0, 0, 1/v^, 0, 0, 0, 0, 0, -l/\/2, 0, 0, 0, 0, 0, 0, 0) 
fs = (0, 0, 0, 1/ V2, 0, 0, 0, 0, 0, 0, 0, 0, -1/v^, 0, 0, 0) 
f4 = (0, 0, 0, 0, 0, 0, 1/V2, 0, 0, -1/^2, 0, 0, 0, 0, 0, 0) 
fs = (0, 0, 0, 0, 0, 0, 0, 1/V2, 0, 0, 0, 0, 0, -1/^2, 0, 0) 
fe = (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -l/A 0) 
ai = (1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)*, 
a2 = (0, 1/^/2, 0, 0, 1/v^, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)*, 
ag = (0, 0, l/\/2, 0, 0, 0, 0, 0, l/\/2, 0, 0, 0, 0, 0, 0, 0)*, 
a4 = (0, 0, 0, l/\/2, 0, 0, 0, 0, 0, 0, 0, 0, l/V^, 0, 0, 0)*, 
as = (0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)*, 
ae = (0, 0, 0, 0, 0, 0, l/\/2, 0, 0, l/V^, 0, 0, 0, 0, 0, 0)*, 
ay = (0, 0, 0, 0, 0, 0, 0, 1/v^, 0, 0, 0, 0, 0, l/\/2, 0, 0)*, 
ag = (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0)*, 
ag = (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, l/v5, 0, 0, 1/^2, 0)*, 
aio = (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1)*, 



Ofi and Oio are 16 x 6 and 16 x 10 null matrices. From 



we have 
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From the singular value decomposition of An we have 



i3i 
B2 



^(0,0,1,0)*, 
^(1,0,0,0)*, 



(l+7b)^l+Dl 



(l+7b)^l+Dl 



(0,D+,0,1)*, 

:(0,D_,0,1)*, 



and 



^1 = -73 (0' 0' 0' 0' 0' 0' 0' 1' 0' 0' 0' 1' 0' 0' 1' 0)*' 
C2 = -75(1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)*, 

(0,0,0,0,0,0,6D+ + i±^,0,0, 



C3 = 

where 



^A_(1+D2) 

1 



(0,0,0,0,0,0,6D_ + i±^,0,0, 



^,6(1 



,0,0,6(1 
.0,0,6(1 



1 + 6 + 662 ± ^1 + 26 + 62 + 2063 + 406^ 



),0,6Z}+ + i±^)*, 
),0,6Z}_ + i±^)*, 
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1 + 6 - 262 ± ^1 + 26 + 62 + 2063 + 406^ 
26(1 + 36) ■ 



Using (j2s|) we have 



Pb 



1 



2(1 + 76) V 1 



/ 1 \ 
10 10 
10 1 
\ 1 / 



;i + 76)(l + D2) ^ 1 
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(30) 



4 Separability of bipartite mixed states 

In this section we discuss some properties related to the Hermitian tensor product decom- 
position that could give rise to some hints to the separability problem of bipartite mixed 
states. From Theorem 3 we can always calculate the tensor product decomposition in terms 
of Hermitian matrices for a given density matrix A, A = Yll=i Bi® Ci. Nevertheless the 
Hermitian matrices Bi and Cj are generally not positive. They are not density matrices 
defined on the subspaces Hi and H2. Hence one can not say that A is separable. 

Let m{A) and M{A) denote the smallest and the largest eigenvalues of a Hermitian ma- 
trix A. We can transform the decomposition into the one given by another set of Hermitian 
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matrices which all have the smallest eigenvalue zero, as follows, 

r r 

A = ^5i®a = ^(Bi-m(Bi)/^ + m(5i)/^)(8(Q-m(Ci)/n + m(Ci)4) 
i=i i=i 

r r 

= Y^{Bi - m{Bi)Im) ® {Ci - m{Ci)In) + Y,MCi){Bi - m{Bi)Im) ® /„ 

1=1 i=l 

r r 

i=l i=l 

r 

= J2(Bi-m(Bi)IJ®(Q-m{Q)Q 

~ r / r \ 

Y,m{Ci){Bi - m(Bi)IJ - m f J]m(Ci)(S, - m(Bi)IJj I„ 

r / r 



(31) 



1=1 



+ 



j=i 



+ 



m 



m 



J2m{B,)cA -J2m{B,)m{a) 



Im ® In 1 



i=l / \i=l / i=l 

where 1^ and /„ stand for m x m and n x n identity matrices. The coefficient of 1^® In-, 
qB,c = m lj2HCi)Bi ) + m I J^^W^i I -J^^WH^i) 



j=i 



1=1 



associated with the decomposition A — Yll=i ® ^ot necessary positive. 

qB,c is decomposition dependent. Associated with another decomposition A = ^[^^ -^i® 
Cj' one would obtain gs',c" 7^ Qb,c- We define the maximum value of qB,c to be the sepa- 
rability indicator of A, S{A) — max{qB,c) for all possible Hermitian decompositions of A. 
With respect to S{A) the associated decomposition is generally of the form 



A = J] 5, Q + /,„ C + S 7„ + S{A)Im ® In, 



(32) 



where Bi > 0, Ci > ,B > 0, C > are positive Hermitian matrices. 

[Theorem 4]. Let A be a Hermitian positive matrix with tensor product decompositions 
of Hermitian matrices like A — Y7i=i Bi® Ci. A is separable if and only if the separability 
indicator S{A) > 0. Moreover S{A) satisfies the following relations: 



S{A) < m{A), 

1 

S(A) > -J2[M{Bi)m{Q) + M{Ci)m{Bi) 

-\m\Bi)\{M{Q) - m(Q)) - |m(C,)|(M(A) - m(5,))], 

S{A) > m{A) - J] M{B,)M{Q). 



(33) 
(34) 

(35) 
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[Proof]. If A is separable, A has a decomposition A = ® the form ((T)), i.e., 

m{Bi) > 0, m{Ci) > 0. We have 

S{A) > qb^c = fn 

> J2 m{Ci)m{Bi) + m{Bi)m{Ci) - m{Bi)m{C{) 

i i i 

= J]m(5i)m(Q) > 0. 

i 

If S{A) > 0, then the associated decomposition (jH^ is already a separable expression of A. 
From the decomposition ()32|) with respect to S{A), we have 

m{A) > J2 HBi ® A) + m{B) + m(C') + S{A) = S{A). 

i 

On the other hand we have 
S{A) > <(B.c>^"("(C,)B,) + 2I"(™(Bi)C,)-^m(B,Ma) 

which is just the formula (jS^)- 

By using the relations M{B + D) < M{B) + M{D), m{B + D) > m{B) + m(L>), 
m(i? + -D) < m{B) + M{D) for any m x m matrices i? and we have 

m{A) = m{Y,iBi®Ci + I,n®C + B®In + S{A)Im®In) 

< m{Im ®C + B®Ln + S{A)I^ (g) In) + M{J2i Bi (S) Ci) 
= m{B)+m{C) + S{A) + M{^,B,®Q) 

< SiA) + j:^MiB,)M{Q). 

Hence formula ()35p follows. □ 

Generally qB,c with respect to our decomposition A = X]i=i Bi®Ci does not equal to the 
separability indicator S{A). Suppose we have another decomposition A = Yll=i B'i®C[. As 
Bi (and Ci) are defined in terms of the singular value decomposition eigenvectors, they are 
linear independent. We can choose linear functionals ipi such that ^i{Bj) = Sij. Applying 
fi^lto both sides of El=iBi ®Ci = EtiB'i ® we get d = EtifiiB'j)q, i.e., 
Ci E < C[, C[., >. Similarly we have Bi G < B[, B[., >. Therefore any other Hermitian 
decomposition A = Yll=i B'i®C'i can be obtained from our decomposition A = X]i=i Bi®Ci 
in terms of the following transformations 

r r 

B'j = EijBi, C'j = ^ Fijd, 

i=l i=l 

13 



as long as the real matrices E = (Eij) and E = (Eij) satisfy the relation EE^ = 1^. 

The inequalities ()33 p - (jH3j) can be served as separability criterion themselves. For in- 
stance, if the minimum eigenvalue of A is zero, then A is entangled if the right hand side of 
is great than zero. 

5 Conclusion and remarks 

We have developed a method of Hermitian tensor product approximation of general complex 
matrices. From which an explicit construction of density matrices on Hi ® H2 in terms of 
the sum of tensor products of Hermitian matrices on Hi and H2 is presented. From this 
construction we have shown that a state is separable if and only if the separability indicator 
is positive. In principle one can always get a Hermitian tensor product decomposition 
of a density matrix by using a basic set of Hermitian matrices. Our approach gives a 
decomposition with minimum terms (the number of the terms depends on the rank of A22), 
similar to the Schmidt decomposition for bipartite pure states. In example pn|l we see that 
the 8x8 density matrix pb has only 4 terms in the tensor product decomposition. 

In j2Zj an entanglement measure called robustness is introduced. For a mixed state p 
and a separable state ps, the robustness of p relative to ps, R{p\\ps), is the minimal s > 
for which the density matrix (p + sps)/{l + s) is separable, i.e. the minimal amount of 
mixing with locally prepared states which washes out all entanglement. In particular, the 
random robustness of p is the one when ps is taken to be the (separable) identity matrix. 
In this case p has the form p = (1 + t)p+ — tim® In/iTin), where pf is separable, p is 
separable if and if the minimum of t is zero. Therefore the separability indicator appeared 
from our matrix decompositions is basically the minus of the random robustness, up to 
a normalization. Another interesting separable approximations of density matrices was 
presented in |28j. This method, so called best separable approximations, was based on 
subtracting projections on product vectors from a given density matrix in such a way that 
the remainder remained positively defined. In stead expressing a density matrix as the 
sum of a separable part and the identity part, this approximation gives rise to a sum of a 
separable part and an entangled part from which no more projections on product vectors 
can be subtracted. 

The results can be generalized to multipartite states. Let's consider a general /-partite 
mixed state pi,2,...,z on space Hi® H2® ... ® Hi. We can first consider pi,2,...,« as a bipartite 
state on space Hi and H2®...®Hi. By using Theorem 2 we can find the tensor decomposition 
Pi,2,-,i — YllLi Bl®Bf^---\ where B] and Bf^- -^ are Hermitian matrices on Hi and H2®...®Hi 
respectively. The matrices can be again decomposed as = Xljli Bfj ® 

Vz, with Bfj and -Bfj being Hermitian matrices on H2 and H^ ® ... ® Hi respectively. In 
doing so at last we have the Hermitian tensor product decomposition of the form, pi,2,...,« = 
'Yl\=iBl ® B^ ® ... ® B\, where are Hermitian matrices on H^. New decompositions 
can be obtained, pi,2,...,z = Eti B'] ® B'] ® ... ® B'\, where B'] = E^'^, Ef^-Sf , A; = 1, /, 
E^ = (E^j) are the real matrices satisfying Yl'j=iEl^jEf^j...Elj = 5i^i^5i^i^...5i^_^ii. 

For any given decompositions, in terms of the protocol (jH!)), one has Pi,2,...,z = X]I=i 
B'^ ... ® B'\ + qldi ® 162 ® ... ® Idi, where Idi is the identity matrix on Hi, {B'], B'^, 
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B'\) are Hermitian matrices such that m{B'\) = 0, or part of them (but not all) are identity 
matrices. The separability indicator S{pi^2,...,i) is the maximal value of the parameter q for all 
possible positive Hermitian tensor product decompositions. If the parameter S{pi^2,...,i) > 0, 
the state pi,2,...,/ is separable, otherwise it is entangled. 
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